arXiv:1507.08572vl [cond-mat.mes-hall] 30Jul2015 


Quantum Confined Stark Effect in Wide Parabolic Quantum Wells 

Sylwia Zieliriska-Raczynska, Gerard Czajkowski, and David Ziemkiewicz 


UTP University of Science and Technology, Bydgoszcz , 
Al. Prof. S. Kaliskiego 7, 85-789 Bydgoszcz, Poland 

July 31, 2015 


Abstract 

We show how to compute the optical functions of Wide Parabolic Quantum Wells (WPQWs) exposed 
to uniform electric F applied in the growth direction, in the excitonic energy region. The effect of the 
coherence between the electron-hole pair and the electromagnetic field of the propagating wave including 
the electron-hole screened Coulomb potential is adopted, and the valence band structure is taken into 
account in the cylindrical approximation. The role of the interaction potential and of the applied electric 
field, which mix the energy states according to different quantum numbers and create symmetry forbidden 
transitions, is stressed. We use the Real Density Matrix Approach (RDMA) and an effective e-h potential, 
which enable to derive analytical expressions for the WPQWs electrooptical functions. Choosing the 
susceptibility, we performed numerical calculations appropriate to a GaAs/GaAlAs WPQWs. We have 
obtained a red shift of the absorption maxima (Quantum Confined Stark Effect), asymmetric upon 
the change of the direction of the applied field (F —¥ -F), parabolic for the ground state and strongly 
dependent on the confinement parameters (the QWs sizes), changes in the oscillator strengths, and new 
peaks related to the states with different parity for electron and hole. 


1 Introduction 

The effects on optical spectra when an external electric field is applied, known in atomic physics as the 
Stark Effect, evolved very rapidly with the invention and development of semiconductor nanostructures. 
The effects of confinement of carriers overlap with the interaction with the field giving rise to the new 
phenomenon known as the Quantum Confined Stark Effect (QCSE). First reported for Quantum Wells by 
Miller et all. m, ®, the QCSE is continuously attracting the interest of the researches. The references 
PQ -[55] are only a small collection of a very large number of papers studying the properties of various 
nanostructures (Quantum Wells, Quantum Dots, Quantum Rods, Superlattices etc.) under electric field. In 
most of this nanostructures the applied electric field causes a red-shift of the positions of the lowest energy 
states, changes in the exciton binding energy, and lowering the oscillator strengths of the resonances. Here 
we consider QCSE in Wide Parabolic Quantum Wells (WPQWs), of thicknesses in the growth direction of 
the order of a few excitonic Bohr radii of the well material (see, for example, [46]-[50] and references therein). 
The optical spectra of WPQWs show a large number of resonances, which are due to the transitions between 
confined states. The Coulomb e-h potential and different confinements for electrons and holes cause mixing 
of the states with different quantum (confinement) numbers. When additionally an electric field is applied, 
states symmetry forbidden appear in the spectra. The behavior of the positions of the resonances is more 
complicated than in the narrow QWs since the lower states show a red-shift, but some higher states show 
a blue-shift or a zig-zag shape, and their oscillator strengths decrease. The case without the electric field 
was discussed in our previous paper [50]. Using the same model, we are able to include the effects of the 
applied field and obtain the solution in analytical form. As an example, we consider a WPQW with GaAs 
as the optically active layer and Gai-^AUAs as the barriers, where the active layer is of the extension of a 
few excitonic Bohr radii, and the constant electric field is applied in the growth direction, which we identify 
with the 2 axis. 
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Although our investigations deal with the theoretical model of WPQW exposed to uniform electric field it 
is believed that such systems are important due to their controllability and potential applications. Employing 
an external electrostatic field to quantum well allows one for steering the optical properties of the system. 
Together with the geometric characteristic of QW the external field is one of the strong modulating factor 
influencing the energy spectrum of charge carriers. Due to controllability of the field the optical properties 
of the nanostructures can be changed on demand. Performing the manipulations of the external interaction 
on WPQW gives one possibility of an effective processing of electrosusceptibility, which may in the future 
be exploited for constructing electrooptical modulators or optoelectronic processors. 

Our paper is organized as follows. In the section 2, we present the assumptions of considered model and 
solve the constitutive equation with effective electron-hole interaction potential and the applied field. Section 
3 is devoted to the details of the applied potential. Next, in section 4, the derived solution of constitutive 
equation is used to obtain the energy levels of the considered GaAs/Gai_ x Al x As wide parabolic quantum 
well. Finally, in section 5, the electrosusceptibility for such nanostructure is calculated and discussed. 


2 The Model 


We will compute the linear optical response of a WPQW to a plain electromagnetic wave 

Ei{z,t) = E i0 exp(i k 0 z - i wt), k 0 = —, (1) 

c 


attaining the boundary surface of the WPQQW active layer located at the plane z = —i/2. The second 
boundary is located at the plane z = i/2. The movement of the carriers in the z direction is determined 
by one-dimensional parabolic potentials, characterized by the oscillator energies hu} e ,huih, respectively. Ad¬ 
ditionally, an external electric held F is applied parallel to the z axis. We adopt the real density matrix 
approach (RDMA) to compute the optical properties (see, for example, 52] - [54] ). In this approach the 
linear optical response will be described by a set of coupled equations: two constitutive equations for the 
coherent amplitudes Y v { 17 , 17 ), v = H,L stands for heavy-hole (H) and light-hole exciton); from them the 
polarization can be obtained and used in Maxwell's Held equations. Having the Held we can determine the 
QW electroptical functions (reflectivity, transmission, and absorption). 

Thus the next steps are the following: We formulate the constitutive equations. The equations will be 
then solved giving the coherent amplitudes Y. From the amplitudes we compute the polarization inside the 
Quantum Well, the electric field of the wave propagating in the QW, and the optical functions. 

The constitutive equation for the coherent amplitude Y in a WPQW and with the applied homogeneous 
electric field F = F k has the form (see, for example, [54] ) 


E q — Turf — i r + 
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Y = M(r)E(R), 


( 2 ) 


where p = yj(x e — Xh) 2 + (y e — Uh ) 2 is the two-dimensional e-h distance, V e h(p, z e , Zh) is the electron-hole 
interaction potential, M(r) is the transition dipole density, which form we have assumed as 


M(r) = M(p, s, $) = ^-8{z)8 (p - p 0 ) 

27rp 0 


(3) 


z = z e — Zh being the relative coordinate in the z direction, po is the coherence radius (the physical meaning 
was explained, for example, in [55], [53]), R jest is the excitonic center-of-mass coordinate, E(R) is the 
electric field vector of the wave propagating in the QW; and p p , P|| are the momentum operators for the 
excitonic relative- and center-of-mass motion in the QW plane. In the consideration of narrow QWs (with 
extension less than one excitonic Bohr radius and arbitrary confinement shape) the following approximation 
was often used. The movement in the 7 - direction was decoupled from the movement in the xy plane, and 
the electron-hole interaction was assumed in the 2 -dimensional form 


V eh = — 
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47re 0 eh p 


(4) 
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with the QW material dielectric constant £&. This approximation enabled to obtain analytical solutions 
for the electron and the hole wave functions, and thus the calculation of the optical properties (see for 
example. |54|b Such metod cannot be used in the considered case of wide QWs (the extension of several 
excitonic Bohr radii) since the e-h interaction retains its 3-dimensional character. As was pointed in Refs. 
32]; ESS) the direct numerical solution of the constitutive equation fl2J is, at the moment, not available 
because of lack of the appropriate orthonormal basis to use in order to decrease the dimension of the 6- 
dimensional configuration space, |51j . Therefore we use the following 3-dimensional form of the interaction 
potential 


V eh 


—S exp 



\2 2 
z h ) -wp 


(5) 


with parameters v, w appropriate for a given nanostructure, which enables to perform analytical calculations 
and reproduces the basic properties of the exciton mm)- 

In the following we assume that the propagating wave is linearly polarized in the x direction, and that 
the vector M has a non-vanishing component in the same direction. We find in the equation ([2]) Hamilton 
operators for the one-dimensional harmonic oscillator 
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with an analogous expression for Hh . Therefore we look for a solutions Y in terms of the eigenfunctions and 
eigenenergies of the operators H ei Hh'- 
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with analogous expressions for the hole, where 
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Hj{x) are Hermite polynomials and Nj normalization constants. Taking into account the valence band 
structure and heavy- and light holes excitons, we obtain 
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where 

and F\h is the so-called ionization field 
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For the GaAs data we have Fjh = 2.318 kV/cm and Fjl = 3.286 kV/cm. Analogously, we have 


1 


m e 


F 


2«e \P\\hJ Fih 

The expression for the Stark shift in the equations (j8]) takes now the form 
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For the light-hole excitons we obtain quite analogous expressions, substituting the respective parameters. 
Using the above functions, we seek the solution for Y in the form 


N N 

^ (P)£e;£h) = 'y ' 1pej{^,e)' l l , nh{^h)Yj n (p) = ^ ) lj n (p) |ejhn), 
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Substituting (fl6l) into the eq. © we obtain equations for the functions Yj n 
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(17) 


We assume the so-called long-wave approximation and consider E(R) in the equation OH) as a constant 
quantity. 

Using the dipole density ©, the model potential ©, and neglecting the center-of-mass in plane motion, 
we put the constitutive equation © into the form 


where 
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Vrsnj = S(erhs exp -v (z e - Zh) enhj). 


(19) 


When the electric field is absent, only states of the same parity will give non vanishing elements (r|s). For the 
field F ^ 0, due to the displacement between the electron and hole confinement eigenfunctions, all possible 
combinations, for example |0e0h), |0e2 h), |le3h), but also |le0h) etc. have to be taken into account. 


3 Calculation of the constitutive equation coefficients 

We have to compute the elements (r\5(z e — Zh)s) and the potential matrix elements 
(er\S(z e - z h )hs) = J J d^d^Vw^tW^O^e - z h ) 

= « / / W-lw - «<)*,>« - a M z. - z h ) 
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Table 1 : Band parameter values for GaAs, AlAs, and Gao^Alo^As, AlAs data from j^H], for Gao.7Alo.3As by linear 
interpolation. Energies in meV, masses in free electron mass mg, 71,72 are Luttinger parameters 


Parameter GaAs AlAs Gao.7Alo.3As 

E, 


m e 

7 1 

72 

*nh\\H 

m h \\L 

M|| H 
Mil L 
Wl'hzH 
rrihzL 
Rh 

R*l 

R*e 

a* H 

a L 

a* e 

tb 


1519.2 

0.0665 

6.85 

2.1 

0.112 

0.210 

0.042 

0.05 

0.38 

0.09 

3.64 

4.3 

5.76 

15.78 

13.265 

9.97 

12.53 


3130 

0.124 

3.218 

0.628 

0.26 

0.386 


0.51 

0.22 

13.32 

19.35 

7.03 

4.84 

11.16 


2002 

0.084 


0.39 

0.13 


12.12 


To perform this calculations, we use the transformation 
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Now the matrix elements (flT)l) can be expressed by integrals 
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Finally 


see, for example, [55]. In particular 
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4 The solution of the constitutive equation 

As we noticed in Ref. m, in order to account the lowest exciton state we take the single function 

Mp,^) = ^=e~ Ap2/2 , 

being the normalized eigenfunction of the Schrodinger equation with the Hamiltonian 

d 2 _ l_d_ 
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where ^0000 = Voooo/A*, zu = a* 2 w. Note that now Voooo depends on the field strength F and so also the 
corresponding eigenvalue £q- We put the function (128[) . with the computed value A, into the expansion (fl 6 l) 
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where now Yj n are constant coefficients. Equation (1181) takes now the form 
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After rescaling the spatial variable p —> p/a* we obtain from (1511) the relation 
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and, in consequence, 


( k rs + e o) Yrs + ^oooot— — Y ra - — p — 'S^'v rsnj Y nj = EM 0 (er\S(z e - z h )hs)i/; 0 (p 0 ) . (34) 
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nj 

We obtained a system of linear algebraic equations for the coefficients Y n j . Having them, we determine the 
amplitude Y (or amplitudes, when accounting the heavy- and light hole excitons H and L. Given the ampli¬ 
tude, we compute the polarization inside the WPQW. For the further calculations we introduce dimensionless 
quantities y rs 

2Mn , N 

- °—Y rs =y rs -E 35 

e 0 e b na* 

and arrived to the formula 

(k 2 rs + e 0 ) y rs + ^ooooy—-— y rs - —^— V' v rsnj y nj = ^^(er\S(z e - z h )hs)ip 0 (po), (36) 

A + w A + w L ' lx 

nj 

where we used the relation 2 = 0- - = —pjr, with Alt being the longitudinal-transversal splitting energy, 

(see for example. [54j ). 




(a) (b) 

Figure 1 : The imaginary part of the mean susceptibility form the eq. ( 1411 ) . for two GaAs/Gao.7Alo.3As WPQWs, a) 
of thickness 51 nm, b) of thickness 32,5 nm 


5 Results for GaAs/Gai-^Al^As Parabolic Quantum Well and dis¬ 
cussion 

The calculation of the WPQW electrooptical functions consists of several steps. First, we define the 
confinement energies Tiuj et h- To this end we choose a specific WPQW having in mind the experimental 
results of Miller et al.[56l. They obtained optical spectra for GaAs(Well)/Gao. 7 Alo. 3 As (Barrier) QWs 
of three thicknesses: L = 51 ± 3,5 nm, L = 32,5 ± 3.5 nm, L = 33.6 ± 3.5 nm. We have per¬ 
formed the calculations for the thicknesses L = 51 nm and L = 32,5 nm. The confinement parameters 
were obtained from the lowest energy levels of equivalent rectangular QWs with confinement potentials 
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V e , con f = 410.38 meV, V h . conf = 72.42 meV gH|. 

The QWs energy states were obtained by standard methods (see, for example, Ref. EZ! and [49] for calcu¬ 
lation details), using the band parameters from Tabic [Q In the considered parabolic QW, the lowest energy 
levels are very similar to the case of an infinite, rectangular well, as one can see on the Fig. [2] 



Figure 2: Comparison of the lowest energy levels in an infinite and parabolic QW, calculated from equa¬ 
tion HI), for GaAs(Well)/Gao. 7 Alo. 3 As (Barrier) and thickness L = 32.5nm. 


The values a*, R* are appropriate for electrons and holes for the QW material, and are defined as 


R* = 


me 


2(4ire 0 e b ) 2 h 2 ' 


a = 


h 2 (47re 0 e b ) 


(37) 


The corresponding values, listed in Table [T] were obtained by using in (E3 the appropriate effective masses: 
m e for R*, a*, and ^\\h,l f° r R*hi a *H an d a L i M|| h,l are the in-plane reduced masses for the electron-hole 
pair and for the heavy- and light-hole exciton data. 
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Figure 3: The imaginary part of the mean electrosusceptibility for the GaAs/Gao. 7 Alo. 3 As WPQW of thickness 32.5 
nm, for different energy intervals 


The results for the confinement energy states are displayed in Table [2] From this energies we obtained 
the confinement energies as 

huj e = 2E e0 , hwhH,L = 2E 0 zH ,l- (38) 
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Figure 4: The same as in Fig. [3] for the WPQW of the thickness 51 nm 


Table 2: Confinement parameters for the WPQWs from Ref. [56], dimensions in nm, energies in meV, 


L 

Tiu) e 

TlLOhH 

hu h L 

E e o 

EozH 

EqzL 

a e a* H 

othHa* H 

&hLa*L 

32.5 

81.66 

19.74 

108.4 

40.83 

9.87 

54.2 

4.21 

4.95 

4.76 

51.5 

43.56 

8.46 

34.4 

21.78 

4.23 

17.2 

3.07 

3.08 

2.68 


Using the above parameters and taking into account the lowest 25 confinement states (see Table 0 we 
have solved the eqn. (ESD and obtained the coefficients Yj n from which we have determined the induced 
polarization inside the WPQW by the relation 


N 

P{z) = 2MoMpo) ^2 I ejhn)(z)Y jn , (39) 

j,n=0 

with the notation 

\ejhn)(z) = ip ej (z)ip hn (z). (40) 

Having the polarization, we compute the mean dielectric susceptibility 

N 

X = ne b ip 0 (p 0 )'^2ye{l\l)A/2 (41) 

e=o 


A/2 

where (1|Qa/2 = y / |Q(C)dC, A = -E, and the index £ runs over the 25 states listed in Tabled Then, 

—A/2 

having the mean susceptibility, one can compute, using the appropriate boundary conditions, the optical 
functions (reflectivity, transmission, and absorption). We have computed the electrosusceptibility for two 
WPQWs of thicknesses 32.5 nm and 51 nm. The parameters S = 2.6, zj = 0.154 and v = 0.5 were determined 
with the procedure described in Ref. [50] . Assuming a certain value of the coherence radius po, we have 
determined the lowest excitonic eigenfunction i/jq. Finally, taking a certain value of the damping parameter 
r, we have solved the constitutive equation (1181) . obtaining the coherent amplitudes. From the amplitudes 
we have computed the mean dielectric susceptibility SD. The advantage of the RDMA is that we obtain 
simultaneously the real and the imaginary part of the susceptibility. The results for the the imaginary part 
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Table 3: Confinement states accounted in computation 


\e0h0) |1) 

|e0/il) ->)|2) 

|e0/i2) ->■ |3) 

\e0h3) ->■ |4) 

|e0/i4) |5) 

\elhO) ->• 6) 

|elM) -4 |7) 

\elh2) -)■ 8) 

\elh3) ->■ 9) 

\elh4) ->• |10) 

\e2h0) ->■ |11) 

\e2hl) 112) 

\e2h2) -)■ 113) 

\e2h3) ->• 114) 

\e2M) ->■ 15) 

\e3h0) j 16) 

\e3hl) 17) 

|e3/i2) 18) 

\e3h3) ->■ 19) 

\e3h4) ->■ 20) 

\e4h0) -> 21) 

\e4hl) ->• 22) 

\e4h2) 23) 

\e4h3) -)■ 24) 

\e4h4) ->■ 25) 



Figure 5: The real part of the mean electrosusceptibility for the GaAs/Gao. 7 Alo. 3 As WPQWs of thicknesses 32.5 
nm for different energy intervals and applied field strengths 


of the mean susceptibility of the considered WPQWs are displayed in Fig.l, Fig.2, and Fig.3. In Fig. [L]we 
show the general effect of the applied electric field for two GaAs/Ga 0 . 7 Al 0 . 3 As WPQws. We observe the red 
shift of the resonances, changes in the oscillator strengths, and the occurrence of new resonances due to the 
broken symmetry. The spectra for F = 0 agree well with the experimental results by Miller et. al. f5f)j and 
our previous theoretical results [50]. In Fig. 2 and Fig. 3 we show the obtained spectra in a more detailed 
form, as compared to Fig. 1. 

In all the cases we observe changes in the placement of resonances, and the occurrence of new peaks 
attributed to different symmetries for the electron and the hole confinement functions. For high values of 
the applied electric field the effects are smaller which is due to the decreasing overlap of the electron and the 
hole confinement functions. We also observe that the shape of the spectra changes with the change of the 
direction of the applied field. Using the properties of the RDMA, we also obtained the real part of the mean 
electrosusceptibility, which is displayed in Fig.4 and Fig. 5a. The impact of high electric fields is displayed 
in Fig.5b, where we show the changes in the real part of the electrosusceptibility for the applied fields up 
to 60 kV/cm. Our method allows to determine the energy shift as a function of the applied field. We have 
computed the energy shift for the lowest confinement states. We observe the quadratic Stark shift for the 
lowest state and a more complicated field-dependence for higher states, as is displayed in Fig.6. Finally, we 
show that the energy shift drastically depends on the thickness of the QW (Fig.7), as was also observed for 
narrow QWs (see, for example, (6], |54j and references therein). Since the confinement energy depends on 
the QW thickness as A -2 , we see that for rectagular QW, A E = —CL 4 . The energy shift in the parabolic 
well on the Fig. I6dl also closely follows this relation. 

6 Conclusions 

We have developed a simple mathematical procedure to calculate the clectrooptical functions of Wide 
Parabolic Quantum Wells. Using the Real Density Matrix Approach and a model e-h interaction potential, 
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(a) 


(b) 




(c) (d) 

Figure 6: a), b) The real part of the mean electrosusceptibility for the GaAs/Gao. 7 Alo. 3 As WPQWs of thickness 51 
nm , for different energy intervals and applied field strengths, c) the impact of high electric fields (WPQW of thickness 
32.5 nm) d) The comparison of the Stark energy shift of the two considered WPQWs for the lowest resonance. 


we derived an analytical formula for the WPQW electrosusceptibility, from which another electrooptical 
functions can be obtained. The presented method has been used to investigate the electrooptical functions 
of GaAs/Gai-^Ala-As WPQWs for the case of radiation incidence parallel to the growth direction. We have 
obtained the red shift of the resonances, changes in the oscillator strengths and new peaks related to elec¬ 
tronic transitions forbidden for the case with absent electric field. We also observed the dependence of the 
spectra on the size of the QW and on the direction of the applied field. For the cases where the experimental 
data were available (for example, for WPQWs with F=0), we obtained a good agreement of our theoretical 
results with experiment. We hope that our results may stimulate experiments on Quantum Confined Stark 
Effect in WPQWs. 
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Figure 7: The Stark energy shift for the lowest resonances: (a)for the GaAs/Gao. 7 Alo. 3 As WPQW of thickness 51 
nm and (b) for the GaAs/Gao .7 AI 0 . 3 AS WPQW of thickness 32.5 nm 
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